Critical slowing down and fading away of the piston effect in porous media 
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We investigate the critical speeding up of heat equilibration by the piston effect (PE) in a nearly 
supercritical van der Waals (vdW) fluid confined in a homogeneous porous medium. We perform an 
asymptotic analysis of the averaged linearized mass, momentum and energy equations to describe 
the response of the medium to a boundary heat flux. While nearing the critical point (CP), we 
find two universal crossovers depending on porosity, intrinsic permeability and viscosity. Closer to 
the CP than the first crossover, a pressure gradient appears in the bulk due to viscous effects, the 
PE characteristic time scale stops decreasing and tends to a constant. In infinitly long samples the 
temperature penetration depth is larger than the diffusion one indicating that the PE in porous 
media is not a finite size effect as it is in pure fluids. Closer to the CP, a second cross over appears 
which is characterized by a pressure gradient in the thermal boundary layer (BL). Beyond this 
second crossover, the PE time remains constant, the expansion of the fluid in the BL drops down 
and the PE ultimately fades away. 

PACS numbers: 65.40.De, 47.10. ad, 47.15. G-, 65.20.+W 



Near critical pure fluids are Newtonian, viscous and 
heat conducting fluids. Their diverging compressibility 
allows a fourth energy transport mechanism ^called the 
piston effect (PE), to become prominent llluiyl. The 
fluid contained in thermal boundary layers (BL) formed 
along heated walls expands, generates acoustic compres- 
sions waves which flash back and forth in the finite size 
sample fluid, and induces an adiabatic heating of the bulk 
much faster than it would be by heat diffusion. The char- 
acteristic time of heat equilibration tends to zero to as 
T 3/2 f or a vdw- fluid, where r is the dimensionless dis- 
tance to the critical point (CP) defined as (T — Tc) /Tc, 
T being the temperature of the fluid and Tc its criti- 
cal temperature. This remarkable finite size mechanism 
contrasts with the critical slowing down of heat diffu- 
sion. The homogeneous character of both temperature 
and pressure (except in the thermal BL for temperature) 
are typical features of the PE. Viscous effects are always 
negligible in bulk supercritical fluid although accounting 
for a strong divergence of the bulk viscosity 0] may lead 
to small scale thermal boundary layer gradients. 

In this letter, we address the ID problem of heat trans- 
port in a near critical fluid contained in a porous medium 
. We consider the gravity vector oriented perpendicularly 
to the planar boundaries of the thin sample (1-D approx- 
imation) or that gravity is equibrated by inertia forces as 
in space laboratories. We basically find that pressure 
gradients which build up progressively while nearing the 
CP (first in the bulk and then in the thermal BL) consid- 
erably weaken and slow down the PE which still exists 



in infinite media but ultimately vanishes. The config- 
uration under study is that of two parallel and infinite 
planes distant by L limiting a homogeneous rigid porous 
medium saturated by a near critical fluid initially at rest, 
at thermodynamic equilibrium, at its critical density p c 
and set just above its critical temperature so that r((l. 
This defines the reference conditions. At t' ^ 0, the left 
boundary is heated with a flux 4>(t') while the right one 
remains thermally insulated. In the rest of this paper, 
subscripts c and are used to denote quantities in the 
fluid phase at the CP and in the reference conditions re- 
spectively. 

As supercritical fluids can be described as compress- 
ible, viscous, heat conducting and newtonian van der 
Waals fluids, we consider the averaged transport equa- 
tions over an elementary volume containing a large num- 
ber of pores to be identical to that obtained for an ideal 
gas in a porous medium. The masss conservation equa- 
tion is given in The momentum balance is expressed 
by the unsteady form of the Darcy equation [fj which 
needs to be valid that the pressure gradient at the pore 
scale is small compared to that at the macroscopic scale. 
The energy balance which acounts for the work of pres- 
sure forces is similar to that of pure fluids, since the one 
temperature model as derived in Q leads to a solid con- 
tribution which only modifies the heat capacity and the 
thermal conductivity of the pure fluid. 

At time t' = 0, heat deposited at x' = diffuses into 
the medium and we derive the asymptotic analysis of the 
coupled equations. To this end, we divide the domain 
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into two sub systems: a thin diffusion BL of thichness S, 
5((l when referred to the sample length, and the bulk. 
The characteristic heat flux at the wall is such that the 
averaged density, velocity, pressure and temperature of 
respective first orders of magnitude fj p ,fj u ,fjp and t\t in 
the BL undergo small perturbations with respect to the 
reference conditions. If we denote 9 = qt, <r((l, 6 is re- 
ferred to a time scale which is <; _1 times longer than 
the sample acoustic characteristic time which defines the 
non dimensional time t as indicated below. The non- 
dimensional, rescaled, governing equations in the BL are: 

Sf)p Pe + S^ 1 f] u u z = (1) 
qfj u u e + S-^-^pPz = - (4/3) (s/K) fj a u (2) 

(pC v )* qfjrfe = -|<5 _1 ry u M 2 

fjpP= {\)mf +\t%~ P (4) 

-5- l T-°- 5 frrf z = ipW(9) at z = (5) 

p=u=f=P=0 at 9 = (6) 

u(z,ff) = at z = (7) 

In these equations, dimensionless temperature, density, 
pressure and velocity are defined from the corresponding 
dimensional quantities denoted with a prime as follows: 

. T(x,t) = ^f 1 = l + T + f(z,0)fj T ; 

• p(x,t) = e^fl = i + p(z,0)v P ; 

. P( x ,t) = ^p =P + f,pP(z,9) 

where r, 70, cq are the ideal gas constant per unit mass, 
the ratio of specific heats and the sound velocity in the 
corresponding ideal gas; u (x, t) = u ^ = f\ u u (z, 9). 
The dimensionless boundary layer coordinate is given 

by: z = yj^ S^ 1 = xd^ 1 where x is the dimensionless 
space variable at the sample length scale while the time 
variable t is referred to the acoustic time and is defined 
by t = e is the ratio, far from the critical tempera- 
ture, of the sample acoustic time ^ to the sample viscous 
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diffusion time ^- , vq being the kinematic viscosity of the 
fluid in the reference conditions. 

For usual critical fluids £ is a small parameter of order 
1CP 8 ; larger values correspond to more viscous supercrit- 
ical fluids or to thinner samples. The Darcy number K 
is the intrinsic permeability K' divided by the squared 
characteristic length of the sample and by the porosity 
e p : K = K' ' I ' (e p L 2 ) . K' is proportional to the squared 
characrteristic pore size and the porosity is the ratio of 
the total pore volume to the total volume of the sample. 



The ratio s/K can vary over several orders of mag- 
nitude (typically from 10 -3 to 10 4 ) while s 2 /K is al- 
ways much smaller than unity. Equation 10 is the in- 
ternal energy balance equation which shows the work 
of pressure forces (— (3/2) S~ 1 fj u u z ) and heat diffusion 

(e f p ro ff -i) J S~ 2 VTT~- 05 f zz ) terms. The near critical 
fluid thermal conductivity, A F , is described by the mean 
field theory 11 i.e. A F = ^ = ^ + t-° 5 = t" - 5 

where A' b is its background value and A' F its critical am- 
plitude. 

In we have neglected the solid heat conductivity 
compared to the diverging one of the near-critical fluid 
and the Maxwell approximation [§] was used to derive 
the equivalent heat conductivity of the porous medium 
to give C^S^~) A F = (^-^zf-^j t~ - 5 . Pro is an equivalent 
Prandtl number defined by Pro = ^fp^- where Dto is a 

2e A' F 

reference thermal diffusivity defined by Dto = V- 

o Ep PcC pQ 

where c' p0 = —Jzi r is the specific heat at constant pres- 
sure of the ideal gas. 

One should note that the mean field theory yields a 
t _1 divergence of the heat capacity at constant pressure 
and thus a heat diffusivity which goes to zero as r 5 . In 
the left hand side of ©, the term {pC v )* is the equivalent 
non-dimensional heat capacity at constant volume given 

by (pCv)* = ^=1 + ^i^P s c s , where c„ = ^ = -^r 
is the dimensionless heat capacity at constant volume of 
the near-critical fluid considered as constant because of 
its weak divergence as r -011 when nearing the CP and 
p s c s the dimensionless heat capacity of the solid defined 
by Psc s = 

Equation 10} is the linearized vdW equation of state. 
Equation © is the boundary condition at the left wall 
(z = 0) where < W{9) ^ 1 is the time modulation of 
the flux deposited at this boundary. The dimensionless 
parameter ip in boundary condition (5) is the character- 
istic heat flux defined as <p = ^'" a , x / J where W max is the 
maximum value over time of the deposited flux at the 
boundary. Equation ©is the initial conditions. 

In the bulk region, the governing equations for the per- 
turbed variables p(x, 9), u (x, 9), P (x, 9) and T (x, 9) are 
similar to 111 1411 except the space variable is x, 5 is re- 
placed by 1, the orders of magnitude are respectively, 
?7p, fj u , fjp,fjT and the boundary and initial conditions are: 

f x = at x = 1 (8) 

p = u = f = P = at = (9) 
u(x,6) = Q atx = l (10) 

Matching conditions between the bulk and the BL re- 
gions close the equations, namely 

lim fj F F (z, 9) = lim fjpF (x, 9) , (11) 

z— >oo x— >0 

where F is a physical unknown. 
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Finding an approximate solution of the equations 
first needs determining the orders of magnitude func- 
tions. Boundary condition Q for temperature implies 
that rjTST°- a (p — rj T . The coupling of the thermody- 
namic variables through the equation of state gives 
ffr = fjp = rfjp while the mass conservation equation 
gives <j?7p = 5 fj u . Moreover, as the matching between 
diffusion and the work of pressure forces in the boundary 
layer (diffusion matches with the bigest term in the equa- 
tion) is expressed by fjuS^ 1 = £5~ 2 t~ - 5 t]t, the pre- 
vious equations finally link space and time scales through 
the relation c,5 2 = er 5 and give the order of magnitude 
fl u = £(p of the velocity in the BL The above orders of 
magnitude allow writing the Darcy equation as 



£-1t - 5 P z / 7o 



-4/3 (e/K)u 



K/e((l 



(12) 
(13) 



If condition l|T3l) is fulfilled, which will be supposed 
valid throughout the paper, and if t))t C 2 = (e 2 /K) , 
which will be further commented later in the text, then 
the following hierarchy holds 



,((l((e/K((e 



-1 T 0.5 



(14) 



and the Darcy equation in the BL reduces to P z = 0. 

The BL solution shows that all the variables go to a 
constant at the outer edge of the BL. Considering the 
matching condition ll lll for the velocity, lim fj u u (z, 9) = 

z — >oc 

lim fj u u (x, 9), and looking for a temperature homo- 

geneization time scale (both temperature in the bulk and 
in the BL are of the same order of magnitude) thus gives: 



smaller than the heat diffusion time which goes to infinity 
as t~ - 5 . This critical speeding up of heat equilibration 
contrasts with the critical heat diffusion slowing down. 
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In other words as far as condition t))t c i = (e 2 /if) 3 is 
fulfilled, the PE is not affected by the porous matrix. 
When nearing the CP, r becomes first of the order of 
t c i and a pressure gradient builds up in the bulk phase. 
Therefore, a first crossover to porous effects arises when 
the pressure gradient and the average viscous stress in 
ltl6)l become of the same order of magnitude, i.e. when 



Tel = 



(17) 



When written in dimensional parameters, this crossover 
value becomes 



T'-T c 



K/ c o) 
K 12/3 



4/3 



(18) 



r e i it thus is a universal value, only depending on the 
fluid and solid matrix properties. 

For carbon dioxide v' = 0.66 .10~ 7 m 2 s~ 1 , c' = 
3 .10 2 m.s^ 1 , and for reasonable values of the ratio 
K' jep = 10~ 14 m 2 , the crossover value is at about 0.09K 
away from the critical temperature. One should note that 
for more viscous critical fluids, this cross over value would 
be much farther from the critical temperature. 

The solution of the system of equations when nearing 
the CP is obtained by performing a crossover type solu- 
tion, e.g. by repeating the above analysis when r tends 

to zero while f — t/ (e 2 /K} 2 ^ is kept constant. Doing 
so leads to the following orders of magnitude: 



Vt = Vt = Sipr - 5 and rj u = rj. 



Sip,. 



(15) 



By carrying the latter into the bulk energy balance and 
matching the transient term with the term representing 
the work of pressure forces (diffusion in the bulk being 
negligible), we obtain the BL thickness S = r and thus 
the value of the time scale <j = £t~ 3 / 2 . 

These orders of magnitude allow writing the Darcy's 
Equation in the bulk as: 



r 3/2 



Ug 



7o 3 



(16) 



If condition t))t c i = (e 2 / K) 3 is verified, the following 

hierarchy holds -^§j2 (("^((iHi 1 !" - an< ^ Darcy's equation 
in the bulk phase also reduces to P x = 0. Under these 
conditions, the governing equations are exactly the same 
as those for the pure supercritical fluid [1(1 1 1 ll | and the 
solution in the bulk displays homogeneous temperature, 
pressure and density. The PE time scale <; _1 , i.e. the time 
scale on which temperature in the sample is homogenized 
by thermo-acoustic heating goes to zero as r 3 / 2 . It is thus 



and 



Vt = Vp = Vt = VP 



3/2 



(19) 



as well as to the following time scale parameter <j = Kj e 

2 

and thermal boundary layer thickness S = {e 2 j~K) 3 . 
As beyond the crossover we have per definition 



r 3/2 



((e 2 /K 



which implies 



(20) 



^{{et/r^ 2 



As 9pe = et/r 3 / 2 is the time variable corresponding 
to the PE in the pure fluid, we can conclude that the PE 
time scale in porous media is longer than the one in the 
pure fluid. The effect of the solid matrix is thus to slow 
down the temperature homogenization whose charateris- 
tic time tends to a constant equal to e/K when referred 
to the acoustic time instead of being shorter and shorter 
as in pure fluid. We call the phenomenon the porous 
saturation of the critical speeding up. 
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Combining Ijllof) with the relevant orders of magnitude 
leads to a diffusion equation for density in the BL whose 
solution can be written as follows 



p(e,z) = --A I K T {9- u,z)W{u)du (21) 



where A (f) = f (70/70 - 1) Pr^ 1 f 05 = At °- 5 is the dif- 



fusion coefficient and 

K T (u 



yf^Ajfju \4A{f) 



is a normalized Gaussian diffusion kernel. Going now 
to the bulk equations we obtain a diffusion equation for 
temperature whose solution is: 



T(9,x) 



A 



(f>Cv) 



K T (9-u,x)W(u)du, (23) 



where 



K T (u,x) = 1 + 2 £ (-1)" 

71=1 

x exp (— n 2 7r 2 B(tildeT)u) 
x cos [n7r(l — a:)] 



(24) 



is the kernel associated with temperature and B (f ) = 
(27/I670 (pC v )*) f 3 / 2 is an effective diffusion coefficient. 
The first correction brought by the presence of the porous 
medium appears in the n = 1 term in the series II24H . 

In the limit f((l, Kt becomes the diffusion ker- 
nel Kt (u) = 



exp 



which indi- 



yj -kB(t)u 

cates that there is an effective energy penetration depth 
y B (f ) in an infinite porous medium saturated by a near 
critical fluid, larger than the heat diffusion penetration 
depth on the same time scale. An effective energy trans- 
fer process in an infinite bulk thus exists whose penetra- 
tion depth, yj B (f ) is longer than the diffusion one on 
the same time scale. This energy transfer, faster than 
diffusion in a convection free infinite medium, can be in- 
terpreted as due to the partial reflection of the heating, 
boundary layer generated, acoustic waves on the porous 
matrix. The thermoacoustic origin of this finite size effect 
could be confirmed by looking at the acoustic propaga- 
tion modes as was done for the PE in a pure supercritical 
fluid 

When nearing the CP much closer than r c i, the 
already-mentioned value r ~ t c2 = (e 2 /K) is reached 
for which a pressure gradient appears in the BL (see 
H12fll. Although the small value of the distance to the 
CP of this second crossover for usual supercritical fluids 



(t C 2 ~ 10~ 5 K for carbon dioxyde for L=10 mm) may 
suggest severe limitations to the Navier Stokes approach, 
a new crossover analysis leads to a diffusion equation for 
density in the boundary layer whose solution, 



9 

J K T2 (0 -u, z)W (u) du , 




(25) 



where the associated diffusion coefficient 
A 2 (f 2 ) 



16 .-3/2 
— 7ot 2 
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and diffusion kernel 



K T2 (u) 



xJttA 2 (f 2 )u 



exp 



4A 2 (f 2 ) u 



(26) 



(27) 



indicates that both molecular diffusion and thermo- 
acoustic energy transfers occur at the BL scale. The 
expresion for the expansion velocity which generates the 
heating acoutic field, 



Uoc (9) = li m u (z, 9) = u{9,x = 0) 



z — >oo 

\A 2 (r 2 )f 2 - 1 - 5 W(9) 



1/2 

indicates that the PE fades away as f 2 when nearing the 
critical point which is however beyond any experimental 
capability in normally viscous supercritical fluids. 
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